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A spherical optical model potential (OMP) containing a dispersive term is used to fit the available 
experimental database of cr(6) and ctt for n+^'^Al covering the energy range 0.1- 250 MeV using 
relativistic kinematics and a relativistic extension of the Schroedinger equation. A dispersive OMP 
with parameters that show a smooth energy dependence and energy independent geometry are 
determined from fits to the entire data set. A very good overall agreement between experimental 
data and predictions is achieved up to 150 MeV. Inclusion of nonlocafity effects in the absorptive 
volume potential allows to achieve an excellent agreement up to 250 MeV. 
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I. INTRODUCTION 



During the last fifteen years, a great deal of theoretical attention has been devoted to casting a proper formulation of 
the nuclear mean field at positive and negative energies. A significant contribution to the solution of this problem can 
be considered the work of Mahaux and co-workers on dispersive optical-model analysis . The unified description 
^ I of nuclear mean field in dispersive optical-model is accomplished by using a dispersion relation, which links the real 
■ and absorptive terms of the optical model potential. The dispersive optical model (DOM) provides a natural extension 
of the optical model derived data into the bound state region. In this way a physically self-consistent description of 
the energy dependence of the OMP is obtained and the prediction of single-particle, bound state quantities using the 
I ' same potential at negative energies became possible. Moreover additional constraint imposed by dispersion relations 
. helps to reduce the ambiguities in deriving phenomenological OMP parameters from the experimental data. 
C3 ' Dispersive OMP analysis has been applied to nucleus-nucleus systems , where the energy dependence of the real 
central potential at low energies near the Coulomb barrier has been studied, and contributions of the dispersion terms 
are evaluated. However for a nucleus-nucleus system, the dispersive OMP analysis is limited to the positive energy 
region, because it is not yet clear how to deal with particle clusters bound in a nucleus. Some progress have been 
achieved in applications of the dispersive OMP analysis to the alpha-nucleus scattering, improving our knowledge of 
the alpha cluster effective interaction inside nuclear system |Q . On the other hand a great success has been achieved 
in deriving DOM potentials for nucleon scattering on closed shell nuclei like ^OCa ^"Zr and ^ospb 

[p]-p|, pl]p^ , p^ , ^ , for which experimental information on bound-states is available. Many studies have dealed also with 
neutron scattering on non magic nuclei(39K ^,^^y ii,^^Kr [||,^^Y ||,"3in (26) and ^''^Bi |2^,|2|,||). 

However very few studies have been devoted to DOM potentials for nuclei with A < 30. Only one preliminary DOM 
^ analysis has been reported for '^''^ Al{n, n ) po). There are two publications making a DOM analysis for proton induced 
reactions on aluminium up to 60 MeV [ |30| , pl[ . 

The main purpose of this contribution is to construct a complex mean field felt by neutrons in ^^Al theoretically 
valid from -50 up to 250 MeV energy. There exist two main versions of the dispersion relation approach. In both 
methods, the real and imaginary parts of the mean field are connected by a dispersion relation and, moreover, the 
mean field is required to closely reproduce the experimental value of the Fermi energy Ep. The main difference 
between the two methods is the following: (i) In the "variational moment approach" the parameters of the 

complex mean field are determined by fitting radial moments of phenomenological optical- model potentials, (ii) In the 
" dispersive optical model analysis" (T^-jl^ , the unknown parameters are derived by performing optical- model fits to 
experimental scattering cross sections that need to be available over an energy range as broad as possible. 
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In the present work a variation of the dispersive optical model analysis is applied to the determination of the 
nuclear mean field for the neutron-^^Al system. An Ohio University - Los Alamos collaboration has published an 
extensive survey of neutron- nucleus total cross section measurements up to 600 MeV |^,^. These high precision 
data together with earlier neutron differential scattering data available in the interval 1-26 MeV form the database 
considered at positive energies. Fermi energy value derived from nuclear masses is used to constrain the mean field 
value at negative energies. Therefore the energy variation of the model parameters is reasonably defined over a wide 
range, an extremely important point for a successful dispersive analysis. Remarkable is the fact that our total cross 
section database goes up to the region where surface absorption can be safely neglected. Since the employed database 
extends up to 250 MeV and since the recent ctt data are very accurate, i.e. the uncertainty Actt is about ±1%, we 
use relativistic kinematics and a relativistic equivalent to the Schroedinger equation in all our calculations. 

Other motivation for our work is that aluminium is an important structural material for the accelerator-driven 
systems and its cross sections are often used as references to determine other cross sections [|4j. There exist phe- 
nomenological OMP (in the sense that dispersive relations constrain is not used) describing neutron scattering on 
aluminium up to high incident energy. The LANL high energy evaluation of Chadwick et al pq] employed the OMP 
of Petler et al ^ up to 60 MeV and the Madland global OMP ^ from 60 up to 150 MeV. Lee and coworkers H] 
derived a new phenomenological OMP which described neutron scattering from ^^Al up to 250 MeV incident energy. 
Recently a new global phenomenological parametrization valid from 1 keV to 200 MeV for A> 27 nuclei was proposed 
by Koning and Delaroche ]39| ]. 

Usually in DOM analysis the absorptive potentials are considered symmetric about the Fermi energy Ep and non- 
zero in the energy gap surrounding Ep. However in the contribution from Mahaux and Sartor they pointed out 
that (i) due to nonlocality effects, the absorptive potential will be highly asymmetric(with respect to Ep) and (n) 
there should be an energy gap centered about Ep in which the absorption term drops to zero, at least for energies 
between the first-hole and first-particle state. Recent DOM analysis of neutron scattering on ■^°*Pb and ^°^Bi |2^] 
failed to describe ax data for energies above 40 MeV using asymmetric version of the absorptive potentials for large 
positive and large negative energies. We will present strong evidence to favour asymmetric absorptive potentials for 
proper description of the neutron scattering ap data for energies between 150 and 250 MeV. 

The paper is structured as follows. Section II provides a description of the dispersive optical model formalism, 
the solved wave equation and the forms of the energy and radial dependencies of the real, imaginary and spin-orbit 
potentials. Section III describes the compound nucleus (CN) calculations, the ^^Al(n,n) experimental database, our 
procedure for searching, and the resulting relativistic and non-relativistic spherical DOM potentials for ^''Al(n,n). 
In the same section we compare derived DOM potentials with phenomenological potentials and experimental data. 
Finally Section IV contains our conclusions. 



II. DOM FORMALISM 



A. Optical-model potential and wave equation 

The optical-model analysis was carried out with a semirelativistic generalization of the conventional nonrelativistic 
Schroedinger formulation of the scattering process |4^. Relativistic kinematics was used for the projectile, but it 
was assumed that the target motion in the center-of-mass system could be treated nonrelativistically. A relativistic 
equivalent to the Schroedinger equation was generated by appropriate reduction of the Dirac equation for a massive, 
energetic fermion (mass m and cm. wave number k) moving in a localized central potential V{r) taken as the time-like 
component of a Lorentz four-vector. In the reduced two-body problem with relativistic projectile but nonrelativistic 
target (mass M) the large component of the partial wave function Fi{p) can be shown to satisfy the radial equation 
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Fi{p)^Q (1) 



where p = kr, is the total cm. kinetic energy, I is the orbital angular momentum, and V{p) is the renormalized 
nuclear optical potential 

V{p)^^U{r),^^l + —^ (2) 
i c + 2m 

Equation (|l|) is formally identical to the radial equation for the solution of the non-relativistic Schroedinger equation 
for the analogous scattering problem with a nuclear potential renormalized by a factor 7. This factor becomes 
increasingly important as the projectile kinetic energy increases (see equation (0)) leading to an effective increase of 
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the potential depth. The spin-orbit term in V{r) employed in this analysis is a purely phenomenological one since 
the intrinsic SO term in the Dirac equation is negligible small in the above limits. Equation was used in all 
calculations. In a non-relativistic case we set a factor 7 equal to 1 and non-relativistic kinematics was employed, 
otherwise relativistic kinematics and the factor 7 according to equation were used. 

Our analysis spans an energy range from 0.1 up to 250 MeV. Both direct and statistical processes contribute to 
nucleon-nucleus elastic scattering at these energies. According to our estimation the statistical processes are important 
up to 12 MeV in aluminium. Compound nucleus calculation will be described in the next section. The direct processes, 
increasingly dominant at higher energies, can be described by the optical model. Although ^^Al nucleus is deformed, 
the spherical OMP has been applied successfully . A posteriori analysis of the impact of this approximation 

on the calculated observables will be discussed below. 

The optical model potential may be written as 

t/(r, E)^- [V^E) + iW^r, E)] f„,{r, R,,a,) 
- [VsiE) + iWsir, E)] g„,ir, i?„ a,) 

[VsoiE) + iWsoiE)] X ^g„^ (r, a,„) (f - a) (3) 

where the successive complex-valued terms are the volume central, surface central and spin-orbit potentials. The 
volume shape f„s{'''i J^d, flu) is a standard Woods-Saxon form factor specified by a potential radius and diffuseness 
a^,. The surface(spin-orbit) shape is the first derivative of the Woods-Saxon form specified by a potential radius 
Rs{Rso) and diffuseness 0^(0^0) 

ffwsC^^^'Oj) = -4a~f{r,Ri,ai) (4) 

The reduced radius parameter is introduced as usual by the relation Ri = riA^I'^ . In our formulation of the OMP 
in Eq. (H) the real and imaginary central volume terms share the same geometry parameters r„ and and likewise the 
real and imaginary central surface terms share the same and as- This assumption |^ can be seen as a consequence 
of the dispersive relations, allowing us to reduce the number of geometrical parameters in the OMP. 
For the spin-orbit potential we adopt the parameters obtained by Koning et al p^ , namely: 

V^o{E) = 6.0exp(-0.005i;) MeV 
W,o(i;) = 0.2-0.011^; MeV 

r,o = 1.017 fm, aso = 0.6 fm (5) 

In a dispersion relation treatment, the real central potential strength consists of a term which varies slowly with 
energy, the so called Hartree-Fock (HF) term, Vhf{E), plus a correction term, AV{E), which is calculated using a 
dispersion relation. The depth of the dispersive term of the potential AV{E) can be written in the substracted form 



AV(E) = - / W(E') ( — ^ dE' 



(6) 



With the assumption that W{E) be symmetric respect to the Fermi energy Ep^ Eq.(^) can be expressed in a form 
which is stable under numerical treatment pq], namely: 



/\V(F\ - p ^ r W{E')-W{E) 



where W{E) is the imaginary part of the OMP. The dispersive term AV{E) is divided into two terms AVv{E) and 
AVs{E), which arise through dispersion relations (|^) from the volume Wv{E) and surface Ws{E) imaginary potentials 
respectively. If imaginary potential geometry is energy dependent, then radial dependence of the dispersive correction 
can not be expressed using a Wood-Saxon form factor, i.e AVv{r,E) ^AV{E)f{r,R,a). However, to simplify the 
problem, the OMP geometry parameters used in this work are energy independent. In this case, using the definitions 
of the equation (^), the real volume Vv{E) and surface Vs{E) central part of the DOM potential are given by 

V,{E)^Vhf{E) + AV,{E) 

Vs{E)^AVs{E) (8) 
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It is known that the energy dependence of the depth Vhf{E) is due to the replacement of a microscopic nonlocal HF 
potential by a local equivalent. For a gaussian non-locality Vhf{E) is a linear function of E for large negative E and 
is an exponential for large positive E. Following Mahaux and Sartor B, the energy dependence of the Hartree-Fock 
part of the nuclear mean field is taken as that found by Lipperheide 



Vhf{E) = VaCxpi-a^^iE-Ep)) 



(9) 



where the parameters Vq and a^^ are undetermined constants. This equation can be used to describe HF potential 
in the scattering regime [Q. 

It is useful to represent the variation of surface Ws{E) and volume absorption potential Wy{E) depth with energy 
in functional forms suitable for the dispersive optical model a naly sis. An energy dependence for the imaginary volume 
term has been suggested in studies of nuclear matter theory [[43| : 



Wy{E)^Ay 



[E~Ef) 



{E-EfY + {B,Y 



where and By are undetermined constants. Following Mahaux and Sartor 
dependence for the imaginary-surface term has been suggested by Delaroche et al 



we adopt n 
[ll5| to be: 



(10) 



4. An energy 



W,{E)=As 



[E - EfY 



{E ~ Ef)"^ + (S, 



&xv{-C,\E-Ef\) 



(11) 



where to = 4 and Aa^ Bs and Cs arc undetermined constants. 

According to equations (^0|) and ( |ll| ) the imaginary part of the OMP is assumed to be zero a.t E — Ef and nonzero 
everywhere else. A more realistic parametrization of Wy (E) and Wg (E) forces these terms to be zero in some region 
around the Fermi energy. A physically reasonable energy for defining such a region is the average energy of the 
single-particle states Ep Q. For aluminium we used a value Ep = —5.66 MeV, obtained by averaging the first three 
particle states reported in the microscopical single-particle level calculation by MoUer and Nix Q] . The experimental 
value of the Fermi energy Ep derived from mass differences is equal to -10.392 MeV. 

Therefore a new definition for imaginary part of the OMP can be written as: 



Wy{E) 



An, 





(E-Ep 



[E - EpY + [ByY 



for Ep <E < Ep 
for Ep< E 



(12) 

and likewise for surface absorption. 

for .Bf < £^ < Ep 

WaiE) = { {E- Epr exp(-a|i? - Ep\) for Ep < E 

{E - Ep)"" + (Bs)"' ^ ' ^ 



(13) 



The symmetry condition 



W{2Ef -~E)^ W{E) 



(14) 



is used to define imaginary part of the OMP for energies below the Fermi energy. Equations (^_2|) and (13) are used 
to describe imaginary absorptive potential in this contribution. 



B. High energy behavior of the volume absorption 

The assumption that imaginary potential Wy[E) is symmetric about E' = Ep (according to equation (|lj)) is 
plausible for small values of \E' — Ep], however as was pointed out by Mahaux and Sartor [Q this approximate 
symmetry no longer holds for large values of — Ef\. In fact the influence of the nonlocality of the imaginary part 
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of the microscopic mean field will produce an increase of the empirical imaginarypart W{r, E') at large positive E' 
and approaches zero at large negative E' Following Mahaux and Sartor we assume that the absorption 

strengths are not modified below some fixed energy Ea- They used Ea = 60 MeV, however this value is fairly arbitrary 
Q . Let assume the non-local imaginary potential to be used in the dispersive integral is denoted by Wy {E) , then we 
can write [|j 



and 



W,{E) = W^{E) 



Wy{E) = WAE)+a 



1 



{Ef-E- Eaf 



[Ep-E-EaY+El 



, for E < Ep-Ea 



E 



[Ef + Eg 

2E 



,3/2 



--^{Ep+Ea 



for E> Ep + Ea 



(15) 



(16) 



These functional forms are chosen in such a way that the function and its first derivative are continuous at E' = 
\Ep — Ea\. At large positive energies nucleons sense the "hard core" repulsive region of the nucleon-nucleon interaction 
and Wy{E) diverges like a^fE. Using a model of a dilute Fermi gas hard-sphere the coefficient a can be estimated 
to be equal to 1.65 MeV"'^/^ assuming that the Fermi impulse kp is equal to 1.36 fm~^and the radius of the 
repulsive hard core is equal to 0.4 fm. On the contrary, at large negative energies the volume absorption decreases 
and goes asymptotically to zero. The non-local imaginary absorption potential Wy(E) and the symmetric imaginary 
absorption potential W{E?) are represented by solid and dotted lines respectively in the lower panel of figure 

The asymmetric form of the volume imaginary potential of equations (|l5|) and ( [l^ ) results in a dispersion relation 
that must be calculated directly from Eq.(^) and separates into three additive terms p6) . Therefore, we write the 
dispersive correction in the form 



t\Vy{E) ^ AT4(i;) + AV<(i;) + AIS(£^) 



(17) 



where tsV^iEi) is the dispersive correction due to the symmetric imaginary potential of equation_(|2|) and the terms 
AF<(iJ) and AF>(i?) are the dispersive corrections due to the asymmetric terms of equations ( |l5[ ) and ([l6|), respec- 
tively. The resulting energy dependence of the dispersive integrals A !/„(£') and A !/„(_£) for both non-local imaginary 
absorption potential W^(i?) and symmetric imaginary absorption potential W(ET) is represented by solid and dotted 
lines respectively in the upper panel of figure ^. While the symmetric case features an equal contribution coming 
from negative and positive energies, in the asymmetric case the negative energy contribution to the dispersive integral 
is very different to the positive energy value. The resulting dispersive correction for the asymmetric case starts to 
increase already for energies above 50 MeV, making a significant contribution to the real part of the OMP. 

It should be noted that non-locality corrections (equations (|l5| ) and (16)) can be used either for volume or surface 
imaginary potential; however, Mahaux and Sartor [Q have shown that nonlocality consideration for the surface 
imaginary potential has a very small effect on calculated cross sections. Therefore in this work we followed Ref. 1^ 
and only considered the effects of nonlocality in the volume absorption. 



III. DISPERSIVE OPTICAL MODEL ANALYSIS 



A. DOM software 



Search optical model codes ECIS95 in the external input mode |47y4q| and COH v 2.2 were used for DOM 
analyses using relativistic and non-relativistic kinematics respectively. A modification was introduced into the later 
code to force equality of the real and imaginary surface (volume) geometry parameters Rs,as{Rv,av) during the 
search procedure, as it is implicit in the equation (0). The code does not include the dispersion relations, therefore 
the dispersion integrals (|^) of the symmetric forms (|l^ ) and ( |l3| ) of the imaginary potential were calculated numerically 
using a Gauss quadrature method [50| ], while the asymmetric contribution was calculated analytically (see Eqs. (16-19) 
of Ref. |46) ) . An auxiliary code system was developed to produce proper input data sets for both optical model codes 
and to calculate for each data set(i.e. for each energy) the quantity according to: 



Ccalc{E, 



^E) 



-1 2 



Aaexp(-B,6l,) 



cxp 



(E) 



(18) 
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Here, f^caidE , 6i)[al'^]^{E)] and aexp{E,9i) [o'*°p(£')], are the difFerential(total) cross sections from the optical 
model calculations and experiments for a given laboratory energy E, respectively, and Aacxp{E, 6'i)[A(j*°p(i?)] is the 
experimental uncertainty reported. The N^- is the number of data points for acxpiE,Oi). Our code system allows to 
finetunc the OMP parameters of interest to minimize the total search of the entire data set. 



B. Summary of the experimental databases 

A survey of the experimental data spanning from 0.1 to 250 MeV used in the DOM analyses is presented in this 
section. The ^^Al(n,n) a{9) data were obtained from Towle and Gilboy |^ at 1, 2,3 and 4 MeV; Tanaka et al fs^jat 



4.8,6,7 and 8 MeV; Kinncyand Percy Q at 5.4,6.4,7.5 and 8.6 MeV; Dagge et al |54| at 7.62 MeVjVelkley et al |5£ 
at 9 MeV; Boerker et al at 10.2 MeV; Whisnant et al Q] at 11,14 and 17 MeV; Petler et al ^ at 18,20,22,25 
and 26 MeV; Bratenahl et al Q at 84 MeV; Salmon [||| at 96 MeV and Van Zyl et al fpl at 136 MeV. The 27Al(n,n) 
Ay{e) data were obtained from Dagge et al iQ at 7.62 MeV and Martin and Walter ^oj at 14 and 17 MeV. These 
polarization data were used only for testing spin-orbit interaction. Energy-averaged total cross sections or for ^''Al 
were obtained from Finlay et al ||3^,Q from 5.3 to 250 MeV. Additional energy-averaged ar data were taken from 
Refs. ^-|7| to be used for comparing predictions of the model. We selected measurements containing several points 
in energy, specially all with data above 20 MeV. In critiquing all the available experimental total cross section data, 
the high resolution cross section data of Ref. was found to be inconsistent with the rest of the data set and was 
ignored in our analysis. 



C. Compound- nucleus corrections 

The statistical model of nuclear reaction according to the Hauser-Feshbach theory |Q with width fluctuation 
corrections as modified by Moldauer j75j is used to compute the CN contributions to the elastic channel. When the 
cross section is averaged over many CN resonances the shape elastic differential cross section can be incoherently added 
to the compound elastic contribution to compare with the experimentally observed elastic-scattering cross section. 
For neutron energies larger than 12 MeV, the cornpqund-elastic contribution can be neglected. The CN cross section 



calculation is built-in inside the search code CoH |49|. Three reaction channels are considered in the statistical- model 



calculations of the "^^Al CN decay: (n,n), (n,p) and (n,a). Transmission coefficients for proton and alpha emission in 



the exit channels are calculated from the spherical OMP parameters by Percy et al [[76| and Arthur and Young |77|(a 
modification of Lemos OMP [Q) respectively. The transmission coefficients in the entrance and inelastic channel 
were calculated using the DOM potential of the present work. 

Discrete level information is used to represent the low-lying states and the Gilbert-Cameron level density formulae 
|f79| are used to represent the high-lying continuum of states. Figure ^ shows the cumulative number of levels as a 
function of excitation energy for the residual nuclei of the three reaction channels. The discrete state data are taken 
from the Belgya compilation contained in RIPL . The vertical lines indicate the cut-off energy between the discrete 
states and the continuum. It is well known that a CN calculation is highly sensitive to the level density parameters 
modeling the continuum of the excited states. We used the "constant temperature" formula |^ to estimate the 
total number of excited states available at excitation energy iJ, N{E) = exp{{E — Eq)/T), where T is the "nuclear 
temperature" and Eq is the energy shift. These two parameters are determined by fitting the cumulative number 
of available experimental states up to some cut-off energy. The level density parameters for all three residual nuclei 
involved in CN cross section calculations are listed in Table |. Cumulative number of levels as calculated by the 
"constant temperature" model using these parameters is shown as solid lines in figure |[ 



D. Search procedure 



It is well known that search routine does not always converge on optimum solution specially when we are dealing 
with strongly correlated OMP parameters. In our DOM analysis we performed a global optimization combined 
with a grid search using a fit in a limited energy region using a maximum number of two fitting parameters 
simultaneously. Our search procedure can be divided in four main steps: 

1. Search for imaginary W^™p{E) empirical potential depth using total cross section data between 70 and 150 MeV, 
neglecting real and imaginary surface contribution. This energy range is selected in order to neglect surface 
absorptive potential in the first iteration. Once empirical values Wy™P{E) were obtained a fit of the absorptive 
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volume potential Wy{E) using equation (|2|) is carried out. In this way volume absorption is fixed, as well as 
dispersive volume contribution AVy{E) to the central real potential, which is calculated by integration. The 
empirical values of the real volume potential depth V^'^p{E) combined with the AVv{E) are used to obtain a set 
of empirical points corresponding to V^™^ (E) . A typical set of empirical values derived in the above described 
way can be seen in figure ^, as obtained with the search code COH. Finally the equation is used to obtain 
the Vq and a^p parameters that offer a best fit to the empirical real potential data. In the fitting process the 
strength Vq was constrained for the DOM predicted first-particle and first-hole states to be centered around the 
experimental value of the Fermi energy. 

2. At each energy for which neutron elastic differential cross section and neutron total cross section data are 
available from 1 up to 26 MeV, we have conducted a best fit by searching on volume real Vy™P{E) and sur- 
face imaginary Wg"^P{E) empirical potential depths. In the first iteration the corresponding dispersive surface 
contribution AVs{E) to the central real potential was calculated by integration from the starting OMP param- 
eters. CN contributions and width fluctuations corrections were considered in all calculation for incident energy 
below 12 MeV. Once empirical values were obtained a fit of the absorptive surface-peaked potential 
Ws{E) using equation ( |l^ ) is carried out. Dispersive surface contribution AVs{E) to the central real potential 
is re-evaluated by integration. The empirical values of the real volume potential depth Vy"^P{E) combined 
with the AVy{E) calculated for these energies are used to increase the set of empirical points corresponding to 
V^^{E). The equation is used to refine the fitting of the Vq and uhf parameters, derived in point (1), 
using the whole empirical set of potential values obtained in steps (1) and (2). We iterate over steps (1) and 
(2) until the emp irical potential strengths were consistent with our predefined energy functional (see equations 
(p|,([T^) and (p^) over the whole energy range 

3. After fixing potential strengths, the optimum geometry parameters were searched for, iterating over steps (1) 
and (2) to redefine the potential strengths corresponding to the optimized geometry parameters 

4. Finally a global x^optimization using the whole experimental database was carried out to obtain the minimum 
in the multiparameter surface. 

E. The ^'^Al{n,n) DOM analysis 

We started our analysis by using non-relativistic formulation to fit the experimental data. Initial values for geo- 
metrical parameters were provided by the energy independent geometry deduced by Whisnant et al ]4l| ] and used by 
Petler et al for phenomenological analysis of the data up to 26 MeV. They found r„ = 1.18 fm, — 0.64 fm and 
rs — 1.26 fm, as — 0.58 fm. Because the general form of the energy dependence of the imaginary potential used in 
the present model is similar to the '^'^Al{n, n) phenomenological OMP of Lee et al ||3^, we used their volume real and 
imaginary potential parameters as a starting point for our analysis. We were using symmetric imaginary absorptive 
potentials according to equations ( p^ ) and (p^), therefore we adjusted 7 parameters, namely (Vq which define 

the smooth energy dependence of the real volume potential and {Ay, By ) and {Ag, Bg^C's) defining the volume and 
surface absorptive potential respectively. After proper values were obtained by this global minimization the energy 
independent geometry parameters were also optimized. The derived non-relativistic DOM potential parameters are 
listed in Table 0. 

The final ctt DOM fits using non-relativistic potential are compared to ^^A/(n, n) data in FigJ|. It should be 
stressed that experimental total cross section data(except the grayed one) shown in this figure were not used in the 
DOM parameter search. We can observe that the experimental total cross section at energies above 130 MeV was 
always underestimated by our nonrelativistic calculations. We can not change the real volume potential depth (or the 
so called Hartree-Fock potential) without spoiling the fits to the differential cross section. One solution could be to 
consider an increase of the radius of real part of the OMP. However this approach would obscure our treatment with 
energy independent geometry. Furthermore, it is theoretically obvious, that relativistic effects and non-locality should 
show up at this energy regime. Therefore we decided to carry out a full relativistic treatment, including non-local 
contribution to the absorptive potential, which will be reflected on the dispersive contribution to the real potential. 

The starting point in this second stage was the non-relativistic DOM potential. We took into account the non-local 
contribution to the volume absorptive potential according to equations (|l^) and (p^. Only one additional param- 
eter was included, namely the energy Ea above which the non-local behavior of the volume absorptive potential is 
considered. In this later minimization the total cross section data up to 250 MeV were included into the experi- 
mental database. All potential parameters changed because of the sizeable contribution of the non-local absorption 
for energies above 40 MeV as can be seen from figure pi In the same figure the total cross section calculated with the 
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non-relativistic DOM potential is shown for comparison. It is interesting to remark that relativistic correction alone 
is clearly not enough for the correct description of the total cross section from 130 up to 250 MeV. The final set of 
parameters of our dispersive relativistic optical model potential is summarized in Table III. 



F. Comparison with the experimental cross section in the energy domain 0.1<E<250 MeV 

We now compare the experimental cross sections with those calculated from our DOM potentials. The geometrical 
parameters of the model and the strengths of the various components are specified in Tables |l| and III . The dispersion 



relations fully determine the real part of the dispersive contribution once the imaginary part of the mean field is 
specified. 

The a{9) relativistic DOM fits are compared to '^'^Al{n,n) data in figure H. In general, the fits to <j{9) are of high 
quality. Very good agreement between experimental data and calculations is observed in the energy region below 12 
MeV, where CN contribution is important. The highest deviation is observed for energies 25-26 MeV located near the 
diffraction maximum. In this energy region a difficulty was encountered during the fit process evidenced by the fact 
that a common set of surface absorptive potential parameters giving acceptable fits to each type of data(differential 
and total cross section) could not be found. The fits to a{9) indicate smaller values of the imaginary surface potential 
depth As parameter while fits to total cross section point to a values larger by about 2 MeV. Experimental (t{9) data 
for energies higher than 26 MeV were not included in the minimization procedure, but our relativistic DOM potential 
displays an excellent agreement with these data. 

The (Tt relativistic DOM fit is compared with the total cross section data and with calculations using phenomeno- 
logical potentials in figures |^ and ||. It should be stressed that only the Finlay et al | |3^^^ experimental total cross 
section data, shown as grey circles was used in the DOM parameter search. In figure the total cross section fit 
is in excellent agreement with the experimental data in the whole energy range from 10 to 250 MeV. The only phe- 
nomenological potential which gives a comparable agreement with experimental data up to 200 MeV is the one by 
Koning and Delaroche being slightly larger than data in the region of the cross section maximum. Madland OMP 
overestimates the experimental cross section by almost 20 % above 100 MeV. The ax relativistic DOM fit is compared 
to the high resolution total cross section data measured by Ohkubo |Q and Rohr et al |Q in figure ||. The total cross 
section fit using relativistic DOM potential is in good agreement with the averaged experimental data in the whole 
energy range from 0.1 up to 10 MeV and practically equal to the cross section derived from phenomenological OMP 
by Petler et al |^^. The total cross section calculated by the phenomenological potential of Koning and Delaroche 
|p9| is smaller than the one calculated by the relativistic DOM potential of the present work in the whole energy 
range from 0.1 up to 10 MeV, but the shape remains quite similar for all compared total cross section calculations. 
Calculation using relativistic DOM potential, including reorientation effects by considering Al as a deformed nucleus 
{13 = 0.4) with ground state spin equal 2.5 is shown as dashed line in the figure ^. This calculation was carried out 
without readjusting any potential parameter to see the effect of deformation on the total cross section. Maximum 
energy in this calculation was equal to the energy of the first excited level to avoid complexity linked to the coupled 
channel approach. We can see that reorientation effects lead to the reduction of the calculated cross section by 10% 
from 0.1 up to 0.8 MeV. The small differences between the solid and dashed curves are a measure of the error incurred 
by the neglect of reorientation effects and nuclear deformation. 

Figure ^ shows the comparison between the experimental analyzing power and differential cross section at 7.5-7.6 
MeV and the predictions of our relativistic DOM. The agreement is good in view of the fact that these data were not 
used in our fitting procedure. It can be seen that CN contribution is still quite important at this energy. Polarization 
measurements at 11 and 14 MeV are compared with DOM calculations in figure |l^. 

Average volume integral for the real part of the optical potential was determined for the relativistic DOM potential 
as well as for the available phenomenological potentials and is shown in figure [l[ In the same figure the "Hartree- 
Fock" , volume and surface dispersive contributions are shown. The biggest difference between our DOM potential 
and the phenomenological ones are located below 50 MeV, where surface dispersive contribution reach minimum and 
then changes the sign, becoming positive. This pure dispersive effect can not be simulated by any variation of the 
phenomenological OMP parameters. It is interesting that real volume integral above 200 MeV is dominated by the 
dispersive volume contribution as a result of the non-locality. 

Average volume integral for the imaginary part of the optical potential was also calculated. In this case they 
are big differences between phenomenological potentials and DOM results as can be seen in figure |l2|. Low energy 
behaviour is different as was the case for the real volume integral, because the dominance of the dispersive contribution. 
However high energy region is also quite different. DOM integral increases with energy as a result of the non-locality 
contribution to the volume absorptive potential. The only phenomenological potential showing similar behaviour is 
the Madland OMP Wf\- His imaginary volume integral goes parallel to the integral calculated using relativistic DOM 
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potential (not considering a discontinuity caused by two different functional forms employed for reduced radius by 
Madland; one below 140 MeV, the second above this value). There is a clear connection between this increase of the 
imaginary volume integral and the saturation of the reaction cross section at energies above 125 MeV as shown in 



figure 13. This behaviour is consistent with the semiclassical estimation of the reaction cross section. The relativistic 
DOM potential reaction cross section reaches a near constant value of 0.3 barn. The asymptotical estimate of the 
reaction cross section is ttR^ equivalent to the reduced radius of 1.03 fm. This value compares well with the averaged 
reduced radius of 1.1-1.2 fm used for the imaginary potential geometry of the DOM potential. It is interesting to 
point out that different reaction cross sections will have a direct impact on cross sections available for any statistical 
model calculations. 



IV. CONCLUSIONS 



In this work we have presented a dispersive relativistic spherical optical model analysis of neutron scattering up to 
250 MeV for ^^Al nucleus. The excellent overall agreement obtained between predictions and experimental data would 
not have been possible without including dispersive terms in the calculations and non-locality effects in the volume 
absorptive potential. New high precision scattering measurements for the aluminium above 30 MeV are necessary to 
establish our analysis on firmer grounds and confirm our present results. 
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TABLE I. Constant temperature level density parameters for residual nuclei in n + Al reaction 



Residual Nucleus EcutlMeV] T[MeV] £;o[MeV] 

^^"Ai in 2.071 -0.678 

^^Na 5.2 1.875 -2.046 

^■^Mg 6.0 2.113 -1.2157 



TABLE IL Optical model parameters for non-relativistic dispersive potential for n -|- Al reaction up to 150 MeV. 



Paramotcr(Unit) Value 

Vq (McV) 52^ 

a^p {MeV'^) 0.0071 

A„ (MeV) 12.5 

B„ (MeV) 58.8 

r„ (fm) 1.20 

ay (fm) 0.65 

A, (MeV) 12.6 

Bs (McV) 3.25 

Cs (MeV-'^) 0.0395 

rs (fm) 1.11 

as (fm) 0.64 

Ef (MeV) -10.392 

Ep (MeV) -5.66 



TABLE III. Optical model parameters for relativistic dispersive potential for n +'^'^ Al reaction up to 250 MeV. 



Parameter(Unit) Value 

Vo (McV) 54~ 

{MeV-'-) 0.0087 

A„ (MeV) 7 

By (MeV) 65 

ry (fm) 1.20 

ay (fm) 0.63 

As (MeV) 12.5 

Bs (MeV) 5 

(MeV-^) 0.034 

(fm) 1.11 

as (fm) 0.64 

(MeV) -10.392 

Ep (MeV) -5.66 
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FIG. 1. Dependence upon energy of the dispersive volume contribution of the real central potential of the n+'^'^ Al mean field. 
The dotted curve corresponds to equation (^o|), in which it is assumed that the imaginary part is symmetric about the Fermi 
energy. The thick solid curves correspond to the asymmetric model, considering non-local behaviour of the imaginary volume 
absorption above certain energy Ea following equations ( p^ ) and (p^. Thin dashed line corresponds to the Fermi energy. 

FIG. 2. Cumulative number of levels as a function of the excitation energy for the three residual nuclei considered in the CN 
cross section calculations. The discrete level data are from the RIPL ^] and are well represented by the "constant temperature" 
level density formula of Ref. using parameters from Table |. Cut-off energy is indicated by the vertical dashed line. Above 
the cut-off energy the "constant temperature" level density formula was used 

FIG. 3. Empirical real volume (solid circles) and imaginary volume potential depth (empty circles) of the OMP for n -|-^^ Al 
as determined from individual best fit searches using atot data in the interval 70 < -E < 150 MeV after the first iteration. 
(Upper panel) The solid line for the Hartree-Fock potential is the functional representation defined in equation (^ . The dashed 
line denotes the starting guess values calculated using Lee et al OMP |Q . The crosses represent the empirical values of the 
Hartree-Fock type potential obtained after the dispersive contribution coming from the volume imaginary part of the OMP was 
substracted from the real volume empirical values. (Lower panel) The solid line for the absorptive potential is the functional 
representation defined in equation ( |l^ ) . The dashed line denotes the starting guess values calculated using Lee et al OMP . 

FIG. 4. Energy dependence of the n -|-^^ Al total cross section from 10 up to 150 MeV. The curve has been calculated using 
the non-relativistic (solid line) DOM potential of the present work. Grey circles correspond to Finlay et al [ p2|j3^ ] experimental 
data used in the fitting procedure. The diamonds, crosses and triangles are obtained from the measurements by Tayloret al 
[ p2[ , Measdayet al and Schneider al |^ 

FIG. 5. Relativistic and nonlocality contribution to the total cross section. The total cross section curves have been 
calculated using the relativistic(solid line) and non-relativistic(dotted line) DOM potentials of the present work. The dashed 
line denotes relativistic DOM potential results without non-locality correction. 

FIG. 6. Comparison between the neutron elastic differential cross section experimental data and our DOM calculations(solid 
line). CN contributions has been added to the direct reaction predictions for incident energies up to 12 MeV. The <7{d) data 
were obtained from Towle and Gilboy jsij at 1,2,3 and 4 MeV; Tanaka et al |^ at 4.8,6,7 and 8 MeV; Kinney and Percy |5^ 
at 5.4^.4,7.5 and 8.6 MeV; Dagge et al^ at 7.62 MeV; Velkley et a/ H at 9 MeV; Boerker et al ^ at 10.2 MeV; Whisnant 
et al Q at 11,14 and 17 MeV and Petler et al Q at 18,20,22,25 and 26 MeV. It should be noted that data above 26 MeV was 
not used in the fitting process. Neutron incident energy is quoted above each calculated curve. 

FIG. 7. Energy dependence of the n Al total cross section above 10 MeV. The curves have been calculated using the 
relativistic(solid line) DOM potential of the present work. The dotted, dotted-dashed and dashed lines have been obtained 
from the phenomenological OMP by Madland [E^, Leeet al |3^ and Koning and Delaroche respectively in their range of 
validity. Grey circles correspond to Finlay et | |32| , ^ experimental data used in the fitting procedure. The diamonds, crosses, 
triangles, empty circles and solid squares have been obtained from the measurements by Taylor al |6^, Measdayet al 
Schneider ef al Franzet al jT^] and Jurenet al |^ 

FIG. 8. Low energy dependence of the n Al total cross section from 0.1 up to 10 MeV. The curves have been calculated 
using the relativistic DOM potential without(solid line) and with(dashed line) reorientation effects. The circles, triangles up 
and triangles down have been obtained from the phenomenological OMP by Harperet al |^^, Petleret al and Koning and 
Delaroche ef al The high resolution experimental data have been obtained from the measurements by Ohkubo [ |n| and 

Rohret al Q 

FIG. 9. The CN corrected o-(9) and Ay (6) data (solid line) at En = 7.62 MeV. Experimental data was taken from Dagge et 
al and Kinney and Perey [g3|. Dashed line is denoted the uncorrected for CN contribution polarization and cross section 
data. 
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FIG. 10. The Ay{0) data (solid line) at E„=14: and 17 MeV. Experimental data was taken from Martin and Walter 



FIG. 11. Energy dependence of the volume integrals per nucleon of the IIartree-Fock(dot-dashed line), volume(dotted line) 
and surface(dashed line) dispersive components of the real part of the n Al mean field. The thick solid curve represents 
the sum of all contributions. Nonlocality was considered in the volume imaginary potential. The solid squares, solid triangles 
and empty circles connected by lines have been obtained from the phenomenological OMP by Madland js^, Leeet al ^] and 
Koning and Delaroche pm respectively. 



FIG. 12. Energy dependence of the volume integrals per nucleon of the volume (dotted line) and surface-peaked (dashed 
line) components of the imaginary part of the n -f-^'^ Al mean field. The solid curve represents the sum of all contributions. 
Nonlocality was considered in the volume imaginary potential. The solid squares, solid triangles and empty circles connected 
by lines have been obtained from the phenomenological OMP by Madland Leeet al |Q and Koning and Delaroche ps| ] 
respectively. 



FIG. 13. Energy dependence of the n + Al reaction cross section from 0.1 up to 250 MeV. The thick solid curve has 
been calculated using the relativistic DOM potential of the present work. The solid squares, solid triangles and empty circles 
connected by lines have been obtained from the phenomenological OMP by Madland [Q, Leeet al and Koning and 
Delaroche |39| respectively in their range of validity. 
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